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ABSTRACT 

We report results of Swift observations for the high mass Be/X-ray binary system lA 
1118-615, during an outburst stage in January, 2009 and at a flaring stage in March, 
2009. Using the epoch-folding method, we successfully detected a pulsed period of 
407.69(2) sec in the outburst of January and of 407.26(1) sec after the flare detection 
in March. Wc find that the spectral detection for the source during outburst can be 
described by a blackbody model with a high temperature {kT ^1 — 3 keV) and a 
small radius (i? ~ 1 km), indicating that the emission results from the polar cap of the 
neutron star. On the other hand, the spectra obtained after the outburst can further 
be described by adding an additional component with a lower temperature {kT ~ 
0.1 — 0.2 keV) and a larger emission radius (i? ~ 10 — 500 km), which indicates the 
emission from around the inner region of an accretion disk. We find that the thermal 
emission from the hot spot of the accreting neutron star dominates the radiation 
in outburst; the existence of both this X-ray contribution and the additional soft 
component suggest that the polar cap and the accretion disk emission might co-exist 
after the outburst. Because the two-blackbody signature at the flaring stage is a unique 
feature of lA 1118-615, our spectral results may provide a new insight to interpret 
the X-ray emission for the accreting neutron star. The time separation between the 
three main outbursts of this system is '^ 17 years and it might be related to the 
orbital period. We derive and discuss the associated physical properties by assuming 
the elongated orbit for this specific Be/X-ray transient. 

Key words: X-rays: binaries — accretion, accretion discs — stars: emission- line, Be 
— stars: neutron — pulsars: individual (lA 1118-615) 



1 INTRODUCTION 

The transient X-ray source lA 1118-615 was first discov- 
ered by Ariel V in outburst during an observation (197 4 
December to 1975 January) of Cen X-3 jEvles et al.lll975h. 
A pul sation of 405.3(6) sec was also detected by lives et al.l 
([1975|) during this observation but it was wrongly attributed 
to an orbital period with an assumption that lA 1118- 
615 consists of two massive compact objects. The optical 
coun terpart was identified as the Be star He 3-640/WRA 
793 (| Chevalier fc Ilovaiskvl [l975|) and is classified as an 
09.5 Ill-Ve massive star with a distance of 5 ± 2 kpc 
(jJanot-Pacheco et al.lll98lh . Long-term and short timescale 
variability i n the Balmer emissio ns of lA 1118-615 has bee n 
reported bv iMotch et al.l (| 19881 ') and iPolcaro et al.l (|l993l l. 
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respectively . This source is also the most extreme Be/X-ray 
transient system known i n terms of the s t rengt h of its Hq 
emission (|Coe et al.l 1 1994 ') . IPolcaro et al.l ([1993) a-lso iden- 
tified a number of other lines (Fen, Hei and ii, Sin, Cm, 
Niii and a tentative identification of Cri); this indicates 
that lA 1118-615 has an extremely unusual Be-type pri- 
mary with complex emission and absorption features. 

The orbital period of lA 1118-615 is not well studied. 
Although the X-ray flux of the source is variable, it remains 
in the quiescent state for most of the time. The Be/X-ray 
binary system underwent a flaring state detected by the 
Burst and Transient Source Experiment (BASTE ) on board 
the C ompton Gamma-Ray Observatory (CGRO) (|Coe et al.l 
Il994 ) an d the WATCH ex periment on board the GRANAT 
satellite (jLund et al.lll992l ) in the late January of 1992. More 
recently, lA 1118-615 showed an outburst that triggered the 
Burst Alert Telescope of Swift (Swift/BAT) on 2009 Jan- 
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Table 1. Log of Swift observations for tlic PC and WT modes 



ObsID Epoch (UT) 



Exposures (s) 
(PC/WT) 



Pulsed counts 
(WT) 



Total counts" 
(PC/WT) 



1 


2009 Jan04 : 


03:36 


105/549 


351 


— /3865 


2 


2009 Jan04 : 


10:05 


— /314 


209 


— /2061 


3 


2009 Jan04 : 


13:21 


5729/12818 


18697 


5729/127740 


4 


2009 Jan06 : 


12:07 


2595/1041 


249 


1969/5950 


5 


2009 Jan07 : 


10:18 


462/2746 


3186 


— /35340 


6 


2009 Mar04 


22:16 


75/4270 


1605 


— /9817 


7 


2009 MarOS 


04:47 


— /5221 


4593 


— /24195 


8 


2009 MarOS 


15:51 


— /3275 


2074 


— /13775 


9 


2009 Marl4 


00:39 


— /5627 


1506 


— /7249 


10 


2009 Marl4 


16:24 


— /4109 


971 


— /5177 



Only the photons observed within 0.3-10 keV were counted. 



uary 4 (IMangano et al.M2009f ). Initial X-ray da ta analysis 
was done by Kond (|2009^ and iManganol (|2009l ). lA 1118- 
615 was detected flaring again during February and March of 
2009 by both the Joint European X-ray Telescope (JEM-X) 
and the soft gamma-ray imager (ISGRI) instruments on the 
Internati onal Gamma-Ray A strophysics Laboratory (INTE- 
GRAL) (JLevder et al.ll2009l ). The X-Ray Telescope of Swift 
[Swift/XKT) also made follow-up observations after the flar- 
ing in March. We here report the Swift observations of lA 
1118-615 in January and March of 2009. 
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Figure 1. Folded 0.3-10 keV Swift light curves of lA 1118-615. 
The upper curve shows the pulsation of 407.69(2) s in the January 
outburst with a pulsed luminosity of 5.3 X 10^^ ergs/s. The lower 
one shows the pulsation of 407.26(1) s after the flare detection in 
March with a pulsed luminosity of 1.6 X 10^^ ergs/s. The d istance 
of the source is assumed to be 4 kpc jCoe fc Pavnelll985r) . 



2 OBSERVATIONS AND DATA ANALYSIS 

We analyzed the X-ray data observed by Swift/yiRT . These 
data were mainly observed within two time intervals. From 
2009 January 4, the Be/X-ray transient was monitored by 
Swift/yiKT for 4 days with a total exposure of 26 ks, which 
includes 17 ks in the Windowed Timing (WT) mode and 
9 ks in the Photon Counting (PC) mode during the out- 
burst. INTEGRAL detected flares of lA 1118-615 again 
from 2009 February 3 to March 3 for a total exposure of 
270 ks. The same field was extensively observed by INTE- 
GRAL and Swift over the following days, and our analyses 
focus on the observations made over an 11 day period with 
Swift/XRT from 2009 March 4 with a total exposure of 22.6 
ks (22.5 ks in the WT mode and only 75 s in the PC mode). 
Details of the observations are listed in Table [l] 

All the Swift/yiRT data were processed with standard 
procedures (xrtpipeline v. 0.12.0), filtering, and screening cri- 
teria by using FTOOLS (v. 6. 5) in the HEAsoft package 
(v.6.5.1). We obtained a refined X-ray position of the source 
by the task "xrtcentroid" from the image of ObsID 1 (Ta- 
ble [l]) investigated with the PC mode. The source located 
at (J2000) R.A.=ll''20"'57^8, decL=-61°54'57".9 with an 
uncertainty of 3". 9 (90% confidence l evel) is the only bright 
X-ray object in the field (|Koneil2009l ). and the source in all 
the other PC images follow this position well. We extracted 
the source counts from a circular region within a radius of 
47" (20 pixels) which is consistent with the 92% of the en- 
circle energy function. The background is obtained from a 
source-free neighborhood of our target with the same ex- 
traction size. 



2.1 Timing Analysis 

The Swift/XHT observations listed in Table [T] were per- 
formed with the PC and WT modes. The data of the PC 



mode retain full imaging and spectroscopic resolution but 
the time resolution is limited to 2.5 s. The data of the WT 
mode only cover the central 8' of the field of view and one 
dimensional imaging is preserved with a better time res- 
olution of 1.7 ms. The observations we used for analysis 
were divided into two groups (January and March) to get 
more counts for accurate temporal analysis. We also ap- 
plied solar system barycentric time correction with the task 
"hdaxbary" using JPL DE200 solar system ephemeris. We 
restricted the data in the effective energy range (0.3-10 keV) 
of Swift/yiRT and per formed epoch-folding with the period 
centered at 405.3 sec (|Nagasg|l989l : lives et al.lll975l ) using 
the resolution of 0.1 sec with 200 trials in the close neigh- 
borhood (~ 395-415 sec). After getting the possible signal 
around 407.0 sec, we then checked the period with more 
delicate resolution (0.01 sec) in the range of ±1.0 sec (406- 
408 sec). The light curve was folded with 32 bins and the 
epoch zero (54835.16972 MJD) was defined at the start of 
the good time interval (GTI) from the ObsID 1. The most 
possible period in our detections is 407.69(2) sec at epoch 
54837.01125 MJD with xl = 152.0 for 31 dof and 407.26(1) 
sec at epoch 54899.95535 MJD with xl = 52.6 for 31 dof. 
(The reference epoch was determined by the mid-point of 
the whole time span.) The corresponding random probabili- 
ties of these two periods are close to 0. The uncertainty level 
of the per i od wa s determined from the x^ value using eq. [l] 
by[Le^|l983) 



^ - o-^i(^^ 



1) 



(1) 



where AP : 



T 



is the inferred Fourier spacing and u is the 



degree of freedom. We show the associated pulse profiles in 
Fig. [l] We note that the normalized strength of the inter- 
pulse around phase 0.2 comparing to the main pulse around 
phase 0.5 increased significantly but the third peak around 
phase 0.8 appears to vanish entirely from January to March. 
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Table 2. Best-fit spectral parameters for the non-thermal spectroscopy of lA 1118-615. The radius is measured from the normalization 
factor for a source distance of 4 kpc and the flux is measured in the range 0.5—10 keV. Pp—teat is the chance probability that the 
improvement of the fit, compared to the single blackbody model (shown in Table [Sjl. 





Observed T 


ime 


Jan. 04 


Jan. 04 


Jan. 04 


Jan. 04 


Jan. 06 


Jan. 06 


Jan. 07 




/model 




03:36(WT) 


10:05(WT) 


13:21(PC) 


13:21(WT) 


12:07(PC) 


12:07(WT) 


10:18(WT) 




Nh (10"" cm-2) 




2 40+0-36 


1 20+0-38 

1-^^-0.31 


1 40+0-22 


1.24" 


9 c-,+0.63 


9 40+0.29 
^-^^^-0.27 


1.08" 


PO 


r 




1.09±0.14 


0.36±-0.17 


c; -3+0. 10 
^-'^•-'-0.09 


0.66" 


n. .0 + 0.20 

U-'i'5-0.18 


0.93±0.11 


0.58" 




Xl/d.o.i. 




1.82/68 


1.61/36 


1.75/100 


3.77/755 


0.86/34 


1.64/104 


2.0.5/514 




Fpo (W-^° erg . 


cm-2 g-i^ 


11.2 


10.1 


9.5 


13.4 


5.9 


8.6 


17.1 




Nh (10"" cm-2) 




, ,q+0.44 
^■■^^-0.39 


-I ,0+0.94 


0.84;o.65 


n cq+O.lO 
^-■^^-0.08 


- 


1 00+0.34 
^•■'''-0.30 


o-42l»:;i 




r 




6.42tl;S 


6.30+2-;s 


4.79tlf, 


4 ie+0.36 
4-15-0.37 


- 


6-92i^-^^ 


5-i8iS:g! 


PO 


fcTpc (keV) 




i-65lS:J^ 


2-02«-it 


9 10+0-1-1 
^-lU-0.13 


, 04+0.03 
1-^4-0.02 


- 


1 ai+0.09 
^•'5 1^-0. 10 


2.00±0.04 


+ 


Rpc (km) 




1.36lO;J| 


0.99«-.i8 


91+0-10 


1.21±0.02 


- 


1 nc;+0.10 

l.UO_Q Qg 


1.31±-0.04 


BB 


Xi/d.o.i. 




1.06/66 


0.77/34 


1.25/98 


1.27/753 


- 


1.18/102 


1.12/512 




FPO+BB (10-10 


erg cm-2 g-l) 


14.2 


14.0 


8.5 


11.2 


- 


9.6 


14.1 




Pp-teat 




0.026 


0.018 


0.28 


1.0x10-10 


- 


0.065 


1.3x10-5 




Observed T 


ime 


Mar. 04 


Mar. 08 


Mar. 08 


Mar. 14 


Mar. 14 








/model 




22:16(WT) 


04:47(WT) 


15:51(WT) 


00:39(WT) 


16:24(WT) 








Nh (1022 cm-2) 




3.57" 


2.25" 


2.78° 


2.63" 


2.85" 






PO 


r 




1.24" 


0.80" 


0.99" 


1.35" 


1.45" 






xi/d.o.i. 




2.61/163 


3.14/348 


2.67/221 


2.61/123 


2.89/90 








Fpo (10-10 erg . 


cm-2 g-i^ 


4.0 


7.0 


6.6 


1.9 


1.9 








Nh (1022 cm-2) 




n 71+0.41 


, OI-+0.23 

1.00_Q 2Q 


, -I < +0.28 


1 or+0.19 

l-'''^-0.18 


9 17+0.51 
^-l'-0.45 








r 




oQ+0.92 
''■■'•'-0.90 


7.77±0.53 


r,,+0.73 


9.50±2.12 


8.11±1.15 






PO 


fcTpc (keV) 




1.58± 0.06 


1 01 +0.04 
l-°^-0.05 


-I 79 + 0.05 


, ^7+0.06 


1 09+0. 08 
^-''^-0.07 






+ 


RpC (km) 




n 07+0. 07 


0.98±0.04 


1 01 +0-06 
1-^1-0. OS 


n cc+0.04 


0.80±0.10 






BB 


x2/d.o.f. 




1.12/161 


1.08/346 


0.95/219 


1.02/121 


1.54/88 








Fpo+BB (10-10 


erg cm-2 g-l) 


77.5 


40.3 


61.5 


11.7 


25.1 








Pp-teat 




1.2x10-12 


6.7x10-33 


1.1x10-20 


1.3x10-11 


3.4x10-5 







The xi °f the fit is larger than 2.0 and XSPEC can not show the uncertainty of the parameter. 



2.2 Spectral Analysis 

Since the center of the source in all the PC mode ob- 
servations is "pile-up" that leads to detector saturation, 
the central source photons has already been eliminated in 
the clean image. In order to avoid the pronounced de- 
pression of counts in the center of our target, we there- 
fore directly downloaded the spectral products of PC 
mode observations from the UK Swift Science Data Centre 
(www.swift.ac.uk/user_objects/). The WT mode observa- 
tions were obtained from the High Energy Astrophysics Sci- 
ence Archive Research Center and proceeded to the spectral 
analysis with the standard data reduction. We downloaded 
the response matrices (rmf) from the official Swift website. 
The ancillary files (arf ) are generated by the task "xrtmkarf" 
with the Swift/yiRT CALDB. Because of a strong absorption 
feature aro und 0.5 keV, which is associated with the Ox ; 
gen edge (|Cusumano fc The Xrt Calibration TeamI I2OO1 



i, 



we only fitted our spectra in the 0.5-10 keV band using 
XSPEC 12.5.1 and rebinned the data with a minimum of 
50 counts per bin to ensure x^ statistics. We produced a 
phase-averaged spectrum for each PC and WT observation 
with enough exposures to generate more than 30 degrees of 
freedom in statistics for spectral fit. Among all the investiga- 
tions, we applied a single power-law model to represent the 
main non-thermal emission of the accretion column/ outer 
gap (Table [2] & Table [4]) and a single blackbody model to 
indicate the thermal emission contributed by the hot spot 



of the neutron star (Table |3] & Table |4|. If the xl oi the 
result is larger than 1, an additional component is required 
to improve the fit. An additional blackbody feature was also 
included in the composite model when we take the thermal 
surface emission of the neutron star (Table [2]) or from the 
accretion disk (Table [3]) into the consideration. The best-fit 
spectral parameters of different models are all given with 
the errors at 90% confidence level. 

For all observations, a single power-law model can not 
provide a reasonable fit as demonstrated in Table [2] (ex- 
cept for the observation of January 06). On the other 
hand, a single blackbody model can provide an accept- 
able fit (Table [3] & (a) of Fig. [2| for the data observed 
in January when we took account the systematic uncer- 
tainties caused by calibrations (Fig. [3}. The spectra with 
long exposures seem to have worse statistics than those with 
short exposures. This might be caused by the complicated 
absorptio n features associat ed with th e response matrices 
(Cusum ano fc The X rt Cal ibration Team 2006). A compos- 
ite model with both power-law and blackbody can be fitted 
with the extracted spectra in March; however, the photon 
index (F ~ 7.5 — 9) is too soft compared with other X-ray bi- 
nary pulsars (XBPs) and the inferred luminosity in the range 
0.3-10 keV will exceed the Eddington luminosity. A compos- 
ite model can also fit most of the spectra in January, but 
the associated photon index of the power-law component is 
still too soft and the main contribution is dominated by the 
blackbody component. We therefore ignore the contribution 
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Table 3. Best-fit spectral parameters for the thermal spectroscopy of lA 1118-615. PC and disk individually represent the related 
parameters of polar cap and accretion disk. The radius is measured from the normalization factor for a source distance of 4 kpc and 
the flux is measured in the range 0.5—10 keV. Pp—test is the chance probability that the improvement of the fit, compared to the single 
blackbody model. 



Observed Time 


Jan. 04 


Jan. 04 


Jan. 04 


Jan. 04 


Jan. 06 


Jan. 06 


Jan. 07 


/model 


03:36(WT) 


10:05(WT) 


13:21(PC) 


13:21(WT) 


12:07(PC) 


12:07(WT) 


10:18(WT) 


Nh (10^2 cm-2) 


0.84«:1I 


47+0-21 


o.49lg:w 


n of, + 0.01 


1 10+0-42 
l-ll'-0.38 


0.97^0-15 


0.29±0.03 


fcTpc (keV) 


1.77±0.09 


9 97+0.23 
^•^'-0.19 


2.i9;°o:lJ 


9 oi+O.Ol 
^■01-0.02 


9 f!i;+0.34 
Z.DO_Q 28 


1.89±0.09 


9 nc;+0-04 
^-Oti_o.o3 


BB Hpc (km) 


1.21±0.09 


„ oQ+0.10 
0.»3_o.o9 


,-, oc+0.06 


1 1C + 0.01 
1.10-0.02 


c:9+0.08 


0-98iro? 


, 9r+0.03 

l-^ti_0.02 


Xl/d.o.i. 


1.14/68 


0.92/36 


1.25/100 


1.35/755 


0.81/34 


1.22/104 


1.17/514 


i=^BB (10-10 erg cm-2 s"!) 


7.3 


8.2 


7.8 


10.8 


4.9 


6.4 


13.8 


Nh (10^2 cm-2) 


1 20 + 0.35 
l-^y-0.31 


- 


n 70+0.42 
0-'«-o.30 


r9 + 0.07 


- 


, 99+0.28 
-^-■^^-0.22 


n of,+0.06 
U. •30-0.04 


fcTpc(fcey) 


1 cc+0.10 


- 


2 lO+O'l-l 
^.tU_Q 12 


1 04+0.03 
l-9'l-0.02 


- 


1 09 + 0.09 


2.01±0.04 


Rpc (km) 


1 Qt; + 0.14 
i.ciO_Q j^2 


- 


nqi+0.09 


1.20±0.02 


- 


1 (14+0.09 


1 90+0.04 

-1-^^-0.03 


2BB kT^iskikeV) 


n TA+0-05 


- 


o.isini 


0.16to-02 


- 


12+0-07 
'J-1^-0.03 


,0+0.03 

U.lU_o.o2 


-Rdisk (km) 


106l?f 


- 


22ilt 


24i^ 


- 


,9, +482 
1^1-101 


109lif 


Xj/d.o.f. 


1.03/66 


- 


1.24/98 


1.27/753 


- 


1.17/102 


1.12/512 


i=^2BB (10-1° ergcm-2 s"!) 


9.5 


- 


8.1 


10.9 


- 


7.3 


13.9 


Pp-test 


0.013 


- 


0.24 


4.6x10-11 


- 


0.046 


1.5x10-5 


Observed Time 


Mar. 04 


Mar. 08 


Mar. 08 


Mar. 14 


Mar. 14 






/model 


22:16(WT) 


04:47(WT) 


15:51(WT) 


00:39(WT) 


16:24(WT) 






N„ (1022 cm-2) 


1 .0+0.16 
J-43_o,i5 


+ 0.08 
0.l't)_o.07 


0.96±0.11 


09+0.14 


009+0.18 






gg kTpcikeV) 


1.77±0.06 


2.05±0.05 


1 qzl+0.06 
l-y*-0.05 


1.59±0.06 


, p;9 + 0.06 

1-0^-0.07 






Rpc (km) 


0.70±0.04 


r, 7n+0.02 
U. '9-0.03 


n 01+0.03 
0-O1-0.04 


r7 + 0.03 


o.eotroi 






xl/d.o.L 


1.55/163 


1.65/348 


1.43/221 


1.52/123 


1.90/90 






Fbb (10-1" erg cm-2 g-1) 


2.6 


5.4 


4.8 


1.2 


1.1 






Nh (1022 cm-2) 


2 c;n+0.34 


1.71±0.19 


1 94+0-25 
l-^*-0.23 


1.21±0.17 


2 06+0--15 






kTpcikeV) 


1 c:n + 0.07 
-■-•3^-0.05 


1 09+0.05 


, -0 + 0.06 
-'--'''-0.05 


, p.p,+0.05 


1.33±0.07 






Rpc (km) 


0.85±0.06 


0.96±0.04 


n qq + 0-05 
^■99-0. 06 


0.64±0.04 


70+0.09 
IJ- '9-0.08 






2BB fc^disk (keV) 


0.14±0.02 


0.14±0.01 


n 14+0-01 


O.O9IO-O2 


0.14±0.02 






^disk (km) 


i79;if 


152lf, 


1891119 


48i;-^ 


,9r + 136 






x2/d.o.f. 


1.12/161 


1.12/346 


0.96/219 


1.05/121 


1.52/88 






F2BB (10-10 erg cm-2 ^-l) 


6.6 


8.6 


9.1 


2.9 


3.1 






Pp-tcat 


2.2x10-12 


3.8x10-30 


3.6x10-20 


8.1x10-11 


1.9x10-5 







January 06: 12:07 (WT) 



March 08: 15:51 (WT) 



1 


j_uuJ0 


%.L 




0.5 


^j^m^ 


\ 


t 


0.2 


7- 


h V 






V 


0.1 






U.Ob 






2 


- 




' ' 


ll ll 1 1 ill 


1 T" 




II 



-2 




■ . 1 1 1 1 ll 


1 ,11 




, , 1 




. .1 


- 



0.5 



1 



10 



0.05 




0.02 



Energy (keV) Energy (keV) 

Figure 2. Spectral fits to blackbody models, (a). The left panel shows the fit to the single blackbody model for the WT mode data of 
ObsID 4 and (b).thc right panel shows the fit to the composite blackbody model (BB+BB) for the WT mode data of ObsID 8. The 
corresponding null hypothesis probabilities of these two fits are all larger than 5% and both of the results are attributed to be acceptable. 
The residuals in terms of sigmas with the error bars of size one are also shown. 
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Figure 3. Spectral fits to blackbody models for 3 spectra observed in January and one spectrum observed in March. The January spectra 
are fitted to a single blackbody component; the March one is fitted to two blackbody components. All of these data have relatively long 
exposures than the other observations and the statistical results of these spectral fits are poor. Some strong absorbed features (e.g. 1.5 
keV of Aluminum, 1.8 keV of Silicon et al.) caused by the uncertainty of the calibrations are clearly visible (the Oxygen edge at 0.5 keV 
had already been removed before performing the spectral analysis). There is also the strong instrumental Nickel contamination present 
in the high energy edge which sometimes is not fully subtracted. If all of these effects can be completely removed, the results of these 
fits will have significant improvements. 



from the power-law for all the spectra and only considered 
a single/composite model with the blackbody component (s) 
to explain the spectral behavior. 

Two-blackbody components can improve the spectral 
fit during the outburst. For ObsID 3 and 5 in WT mode, 
an F-test shows that the additional soft component during 
the outburst is significant at > 99%. But for the other 3 
observations in January, an additional soft blackbody com- 
ponent does not improve the fits significantly. In contrary, 
for the observations detected after the flaring in March, a 
two-blackbody model fits significantly better than a single 
blackbody one in Table |3] ((b) of Fig. [2]). However, we note 
the WT spectra in ObsID 1, 2 & 4 can be described by an 
acceptable single blackbody model with much shorter ex- 
posures than the other observations. It is unclear whether 
the single blackbody fits might be due to the larger statis- 
tical uncertainty in shorter exposures. However, the single 
blackbody component fits with long exposures in January 
are still much better than the single blackbody component 
fits of March observations. This indicates the variation of 
the spectral behavior might be real. 



We also note that the fitting parameters of Jan 04 ob- 
servations (ObsID 1 & 2) are very different. This might be 
caused by the fact that the fitting parameter of the hydro- 
gen absorption is very sensitive to the soft band spectrum 
if the the photon counts are low; a very small variation of 
the spectrum at soft band will result in a complete different 
hydrogen absorption, which will then cause a large variation 
in other fitting parameters. For the spectra of Jan 04 obser- 
vations (ObsID 1 & 2), if we fixed the hydrogen absorption 
to be ~ 6.4 X 10^^ cm-^ the ObsID 1& 2 of the WT mode 
observations can both be fitted by a single blackbody model 
with kT = 1.86±0.07 keV, R =1.11±0.05 km with xl = 1-21 
with 69 d.o.f and kT = 2.lQt°o^^i keV, R =0.89 ± 0.07 km 
with xS = 0.95 with 37 d.o.f. The variations of the parame- 
ters between these two spectra with the fixed hydrogen ab- 
sorption are much smaller than those listed in the Table |31 
This indicates that the huge variations of the polar cap tem- 
peratures and radius within the same day may not be real. 

For all the spectra extracted in March, only the second 
spectrum of March 14 (ObsID 10) does not provide an ac- 
ceptable fit with a two-blackbody model; this might be due 



6 L. C.-C. Lin, J. Takata, A. K. H. Kong and C. Y. Hwang 



Table 4. Best-fit spectral parameters for the pulsed speetroscopy 
of lA 1118-615. The radius is measured from the normalization 
factor for a source distance of 4 kpc and the flux is measured in 
the range 0.5—10 keV. 





Observed Time 
/model 


January 
(WT) 


March 
(WT) 


PO 


Nh (10^2 cm-2) 

r 

xl/d.o.f. 

FpO (10-1" erg cm-2 g-^) 


0.64 (fixed) 

0.23±0.07 

1.23/314 

3.6 


1.88 (fixed) 
, n/i+0'12 

1.28/156 
1.2 


BB 


Nh (10^2 cm- 2) 

fcTpc (keV) 

Rpc (km) 

xi/d.o.{. 

Fbb (10-" erg cm-2 g-l) 


0.64 (fixed) 

1 -70 + 0.10 

0.70±0.05 

0.76/314 

2.8 


1.88 (fixed) 

1.16±0.07 
„ 04+0.09 

0.74/156 
0.83 



to the serious Nickel contamination at 7-8 keV that is not 
fully subtracted from the calibration and many absorption 
features. Some of t he features are caused by the instrumental 
effects as shown bv lCusumano fc The Xrt Calibration TeamI 
((2003); however, the true spectra of the source might be 
much more complicated than our models and the instru- 
mental effects probably do not explain all the poor fits to 
the data. 



2.2.1 Pulsed Spectral Analysis 

To further study the emission from the neutron star, we 
also generated the pulsed spectrum of lA 1118-615 with 
the data observed in the WT mode. We divided each obser- 
vation into pulsed and unpulsed emission depending on their 
pulsed phase shown in Fig. [l] The light curve of lA 1118- 
615 was folded with 32 bins and we assign the phase from 
0.375 to 0.90625 (12th - 29th of 32 bins) to be the pulsed 
phase of the 407.69 s period and the phase from 0.125 to 
0.75 (4th - 24th of 32 bins) to be the pulsed phase of the 
407.26 s one. We also used the same criterion to estimate 
all the pulsed photons in 0.3-10 keV for each data listed in 
Table[T] To get enough counts for pulsed spectral fits, all the 
pulsed spectra in January and March were combined into a 
January spectrum and March one respectively. The associ- 
ated ancillary response files were also combined with weight- 
ings depending on the relative exposures (shown in Table [T]). 
The pulsed spectrum was then obtained by subtracting the 
unpulsed emission from the pulsed emission. According to 
Table [S] the average hydrogen absorption in the best fit to 
a single blackbody spectrum obtained at the outburst stage 
in January is ~ 6.4 x 10^^ cm"^ and is ~ 1.9 x 10^^ cm"^ 
to a two blackbody spectrum after the flare detection in 
March. In the fits to the pulsed spectra, we fixed these val- 
ues to associate with the pulsed-average spectra (Table |4]). 
The difference of the absorption can be caused by the pos- 
sibility that the absorber is related to the different amounts 
and/or the viewing angle of the material accreted by the 
pulsar along the line-of-sight. 

The pulsed spectra of January and March can not be 
well fitted by a single power-law, and we can not find any evi- 
dence of the non-thermal emission from each phase-averaged 
spectrum. However, a single blackbody model can be ap- 
plied to the pulsed spectra in January and March at 0.5- 



^^ ergs cm '^ s """ 
^^ ergs cm^^ s^^ in March 



m 



with xl = 0.76 for 314 d.o.f and kTss = 1.16 ± 0.07 keV 
Rbb = 0.84t2g^ km with xl = 0.74 for 156 d.o.f. The in- 
ferred pulsed flux is (2.8 ± 0.6) x 10 
January and (8.3 ± 2.6) x 10 
The decrease of the luminosity might be caused by the cool 
ing of the surface temperature. 



3 DISCUSSION 



10 keV with kTss = 1.78! 



keV, Rbb = 0.70 ± 0.05 km 



According to A riel V (JEvles et al.lll975h . CGRO/BATSF, 
(|Coe et al.lll994l ) and Swift/BAT observations, the time sep- 
aration of ~ 17 years between the detection of each outburst 
for lA 1118-615 suggests that the neutron star does not in- 
teract vio lently with the mas sive companion during most of 
the orbit (|Villada et al.lll999l ) . If the outburst is only caused 
by the approaching of the neutron star to the periastron, the 
long time separations of the outbursts indicate a very flat 
elliptical orbit. 

However, the orbital peri ods for known Be/XB Ps range 
from several days (A 0538 -66; Johnston et al.ll 19801 ) to hun- 
dred of days (X1145-619; IWatson et al.lll98lh and no other 
convincing evidence can connect this long time separation 
to be the orbital period of lA 1118-615. The time separa- 
tion of ~ 1 month between the outburst of 2009 January 
detected by Swift/XRT and the flare of 2009 February de- 
tected by INTEGRAL might also indicate a possible orbital 
period. One major concern is that we would expect to have 
observed much more flares during the interval of 17 years 
or at least more than one flares after the main outburst if 
the orbital period was only around 1 month. In addition, 
each main outburst detected every 17 years is the normal 
outbursts (type I, 10^® - 10^^ erg/s). However, we would ex- 
pect to observe a gian t burst for a long recu rrence interval 
(type II, > 10^'' erg/s: ICaballero et aLfeoOSl ). On the other 
hand, if the orbital period for lA 1118-615 is 17 years, this 
binary system will have the longest orbital period than all 
the other known Be/XBPs shown bv lCorbet et all (|2009l '). 

We would expect to observe some evidences of the pas- 
sage of the neutron star during the periastron in other wave 
bands. ICoe et al.l (|l994l ) found that the equivalent width of 
Ha varied on very short time scales around the outburst 
of 1992; this can be explained as the disruption of the cir- 
cumstellar disk caused by the passage of the neutron star 
during the periastron. This indicates that the assumption of 
the 17-year orbital period is consistent with the H^ equiva- 
lent width variation. Further optical investigations and the 
study of high-energy emission mechanism of this Be/X-ray 
transient may help us to verify the orbital period and solve 
the related mystery. 



3.1 Observational features detected by Swift 

The pulsation in January was measured during the outburst 
while in March it was detected after the flare detection of 
February. The pulsations of lA 1118-615, shown in Fig. [ij 
obviously demonstrate that the pulsar has spun-up during 
the period between the outburst and flare events. The pulse 
periods do not show variation during the different observa- 
tions of January and March respectively. As a result it is 
unclear whether the outburst or flare is responsible for the 
spinning-up of the pulsar. But the change of the pulse profile 
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for lA 1118-615 may relate to the decrease of the pulsed lu- 
minosity; simi lar luminosity depe ndencies were also claimed 
for C en X-3 (JNaease et a"Llll992h. LMC X-4 (iLevine et all 
Il99ll ) and EXO 2030+375 (|Parmar et al.lll989l ). In the evo- 
lution of the pulse profile for EXO 2030-1-375 following the 
first outburst decay, an interpulse appeared and became 
stronger as the luminosity decreased from the maximum and 
this trend continued until the main pulse had become en- 
tirely disappeared. Then the profile became dominated by 
the contribution from the interpulse. The same luminosity 
dependence was also evident for lA 1118-615. The pulse 
profiles of lA 11 18-615 show featur es similar to those of 
EXO 2030-1-375 (jParmar et al.lll989l '): the pulsed luminos- 
ity decreased and the normalized strength of the interpulse 
became significant relative to the main pulse after the fiar- 
ing. A long term survey is needed to further examine the 
dependence of the pulse profile and X-ray luminosity. 

According to TableO a single blackbody model can pro- 
vide an acceptable fit with a temperature of kT > 1 keV and 
a smaU hot spot radius (~ 1 km; assuming a source distance 
of 4 kpc: ICoe fc Pavnell985l ) for most of the data observed in 
the outburst. The thermal emission can be attributed to the 
ionized photons or the collision of the heated gas; however, 
the most plausible origin is f rom the surface (polar cap) 
emiss i on of the neutron star (ILa Palombara fc Mereghettil 



I2OO7I : iLa Palombara et al.ll2009l '). iHickox et all (|2004l ') pre^ 
sented a clear physical picture to show the emitting process 
during the accretion. In this picture, the radius of the ac- 
creting polar cap is ~ 0.1-Rjvs, where Rns is the radius of 
the neutron star, typically as sumed to be 10 km . In contrary 
to the spectra of most XBPs (jHickox et al.ll2004h . we did not 
find any obvious evidence of a main non-thermal contribu- 
tion. Instead of X-rays emitted from the accretion column 
or outer gap, we detected strong emission from a surface 
layer. This X-ray contribution during the outburst can be 
originated from the thermal emission around the compact 
source, which is similar to the polar cap emission of an iso- 
lated neutron star. 

The second thermal feature [R ^ 100 km & kT < 
0.2 keV) rese mbles the X-rays fr om the inner edge of an ac- 
cretion disk (jHickox et al.ll2004l '). This additional soft char- 
acteristic may become an important component when the 
pulsed flux of the source weakens. Such X-ray thermal con- 
tributions are always obser ved in the bright XBPs with lu- 
minosity > 10^^-^^ ergs/s JPaul et al.ll20o3 : iRamsav etld] 
|2002| : iManousakis et al.ll2009h . but the source luminosity of 
lA 1118-615 in 0.5 - 10 keV observed by Swtft/XRT is only 
~ (0.6 — 1.7) X 10^® ergs s~^. Our results might indicate 
that the soft excess contributed by the emission from the 
inner disk can generally be observed for the XBPs after flar- 
ing, not only for the bright targets. If the variation of the 
second thermal component is real, our spectral analyses in 
different time intervals may indicate that the surface (polar 
cap) emission of the neutron star dominates the observed 
X-rays from lA 1118-615 at the beginning of outburst. The 
fact that the thermal emission contributed by the accre- 
tion disk becomes signiflcant at another specific time inter- 
val may represent that a different violent radiation process 
starts in one month and may be also associated with the 
flare detection. Although the average source luminosity in 
different time intervals only changes by a small fraction, the 
cooling of temperature in the similar polar cap area clearly 



exhibit a decay of the surface emission from the neutron 
star. The decrease of the pulsed emission derived from our 
phase-resolved spectrum also supports this viewpoint. 



3.2 Parameters of the binary system 

The parameters of the neutron star and the properties of the 
orbital motion for the binary system, lA 1118-615, are not 
well understood. Furthermore, the possible orbital period of 
~ 17 yr is extremely large compared with 100 days for the 
typical Be/XBPs with long spin period. Therefore, it will 
be worth examining the parameters of the neutron star and 
orbital motion for further discussion on this unique high- 
mass X-ray binary system. 

The increasing X-ray emission at the outburst stage 
in January of 2009 allowed us to find a rotation period of 
P ~ 407.69(2) sec. On the other hand, a rotation period 
of P ~ 409.2(8) sec was reported by Chandra observation 
on August 21 of 2003 in the quiescence stage of the source 
(jRutledge et al.l 120071 ). Assuming that this spin up of the 
neutron star from its quiescence stage is a result of a contin- 
uous accretion of matter onto the stellar surface of the neu- 
tron star, we can estimate the accretion rate and the mag- 
netic field of the neutron star. The Chandra observation in 
2003 detected an X-ray luminosity of Lx ~ 1.8 x 10"^* erg/s 
in the 0.5-10 keV energy bands. Assuming the X-ray emis- 
sion is from the accretion onto the stellar surface, the ac- 
cretio n rate at the surface is estimated as IjCampana et al.l 
Il998l ) 



M, 



-acc 



Rns r 



5 X 10- 



,13 



'r]BcMi 4 RiiLx;ii 



i/s 



(2) 



where r]BC is the bolometric correction, Rns and Afjvg are 
stellar radius of and the mass of the neutron star, respec- 
tively. In addition, Lx,34 is the X-ray luminosity in units of 
lO"^* erg/s, and Re is the stellar radius of the neutron star 
in units of 10^ cm and Mi. 4 is the mass of the neutron star 
in units of the solar mass. 

We assume that the accretion rate is constant in the qui- 
escent stage before the onset of the outburst in 2009 January. 
Combining the estimated accretion rate Mace with conser- 
vation of the angular momentum that 



ISn/St = MaccVGMNsRn 



(3) 



with Rm being the magnetospheric radius where the mag- 
netic pressure is balanced with the pressure of the accretion 
disk, the strength of the magnetic field of the neutron star 
is given by 



B ' 



2 X 10"r,sSMf.^4'fi6^'ix'34 



(4) 



g ■ cm . 

The strength of the magnetic field can be estimated as 
S ~ 2 X 10" G for M1.4 = 1 and 7? = lO'' cm. This 
value is also similar to tho se extremely magnetized n eutron 
stars (known as magnetar, IWoods fc Thompsonll2006l ) if the 
bolometric correction is tjbc ~ 1- On the other hand, a 
typical strength of the magnetic field of a young neutron 
star B ~ 10^^ G is expected if tjbc ~ 5. Accretion onto 
the neutron star surface is permitted when the magneto- 
spheric radius Rm becomes smaller than the corotation ra- 
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dius, Rcor = (GMArsP^/47r^)^/^, where P is the period of 
the neutron star. This condition is satisfied when r]BC ^ 2. 
It will be expected that the outburst occurs as a result 
of the interaction between the equatorial disk of the Be star 
and the neutron star at the periastron passage. Near the 
periastron, the velocity of the neutron star, vns, will be- 
come vns ^ {GMbb/ RbcY {^ 500 km/s), where we used 
stellar mass and radius of Be star of Mse = 17Mq and 
Rbc = IO-R0, respectively (jSigut et al.ll2009l ). This implies 
that the velocity of the neutron star is in general larger than 
the typical velocity for the matter of the equatorial disk of 
the Be star, whose typical initial vel ocity is << 1 km/s 
(|Marlborough et al.|[l997l : |PorteJll99^ '). The velocity of the 
neutron star dominates that of the matter of the equatorial 
disk of the Be star near the periastron. Then the capture of 
the matter from the equatorial disk of the Be star by the 
neutron star will occur at a radius of 



'2GMns/vI ~ 3.8x10^ 



VNS 



Mns 
IAMqJ VlOO km/s 



(5) 

from the neutron star, where Vr ~ vns is the velocity of the 
material of the equatorial disk relative to the neutron star. 
At the outburst stage in January of 2009, the average X- 
ray emission at 0.5-10 keV was measured with a luminosity 
of Lx ~ 1.5 X 10"^^ erg/s, indicating that the accretion rate 
onto the neutron star is 



Mace ~ 6.4 X lO~'°^M0/yr. 





(6) 



This accretion rate for a young neutron star is about 2 
orders of magnitude smaller than typic al mass loss rate 
(M,„ ~ 10"'^ Mm/vr: IWaters et al.lll988l 'l of an equatorial 
disk around a Be star. We can also estimate the accreting 
density of the surrounding medium from Bondi accretion 
{Mace = 'i'^'TaccpVr', where p is the ambient density). Thus 
we obtain: 



P = 



Ma 



bxlQ-'^riBC 



M. 



NS 



I AM, 



Vns 



100 km/s 



■ g/cm' 

(7) 



This is several orders smaller than the typical value {p ~ 
lQ-10 _ i Q-is g/cm^') of t he density at the inner region of 
the disk (jSigut et al.ll2009l ). implying that (1) the neutron 
star is likely interacting with the outer edge of the Be star 
disk or (2) the disk of Be star in lA 1118-615 is relatively 
thin. 

Our spectral analyses show one thermal polar cap com- 
ponent during the outburst stage in January and show 
an additional thermal disk component after the flaring in 
March. This might be explained as a result of the neutron 
star approaching the equatorial disk of the Be star. When 
the neutron star flrst approached the equatorial disk, the ac- 
cretion disk of the neutron star could have accumulated rel- 
atively little material, so the accretion disk would not have 
significant emission and we could only detect the polar cap 
emission of the neutron star. The flare of February might be 
caused by disk instabilities of the accretion disk of the neu- 
tron star. This indicates that the neutron star would have 
accreted signiflcant amount of material from the equatorial 
disk of the Be star during the periastron; this makes the 
thermal disk component visible in the March observations 
after the February flare. 



4 CONCLUSION 

The Be/X-ray transient lA 1118-615 has been observed in 
outburst by Swift starting in January of 2009. Comparing 
with the epochs of two previous outbursts for lA 1118-615 
detected by Ariel V and CGi?0/BATSE, the time separa- 
tion of ~ 17 years might give a hint as the orbital period. On 
the other hand, the time interval of ~ 1 month between the 
detection of the outburst and the flaring in 2009 might also 
give another indication of the orbital period; especially given 
the fact that the new Be/XBP (e.g. SAX J2103. 5-^4545; 
iBavkal et al.ll2000l ) was discovered to have a rapid orbital 
motion with a long spin period. We performed timing anal- 
ysis on the Swift/yiKT archive data and the pulsations at 
407.69(2) sec and 407.26(1) sec were individually detected 
in January and March. lA 1118-615 had obviously spun- up 
due to the accretion from the massive primary between the 
January and March observations. 

We also investigated the X-ray spectral behavior dur- 
ing the outburst in January and after another flare detec- 
tion in March. Almost all the phase-averaged spectra can 
be well-modeled by thermal emission. During the outburst, 
we found that a single blackbody component with a higher 
temperature of ~ 1 keV and a small radius of _R ~ 1 km pro- 
vides a reasonable fit, and it indicates that most of the ob- 
served X-rays are from the polar cap emission of the compact 
star. After the flare detected by INTEGRAL, the following 
Swift/yiKr observations display a significant soft thermal 
excess above the main surface emission from the neutron star 
in their spectral behavior. The soft thermal excess shows a 
lower temperature of ~ 0.1 — 0.2 keV and a larger emitting 
radius of ~ 100 km; this component might be connected to 
the emission from the inner accretion disk of the binary sys- 
tem. In contrary to other XBPs, we do not find a significant 
non-thermal component in the spectrum of lA 1118-615. 
Our results show that the thermal radiation from the hot 
spot of neutron star and from the reprocessing of the dense 
material in the inner disk provides the emission mechanism 
for lA 1118-615 after the outburst. The highest luminosity 
at 0.5-10 keV of this source detected after the outburst is 
only ~ 2.7 X 10'^^ ergs/s. Our results also demonstrate that 
the soft excess contributed by the inner accretion disk may 
not only be observed in bright XBPs (L^, > lo-^^"-^* ergs/s). 

iHickox et al.l (|2004l ') gave a clear physical picture to ex- 
plain the possible origin of the soft excess. Our discoveries 
based on the spectral analysis may provide a hint to exam- 
ine the accreting phenomenon. We also derived the physical 
properties of this accreting system using a simple model. The 
accretion rate and the density of the circumstellar disk de- 
rived with the assumption of the elongated orbit of lA 1118- 
615 are lower than those of the other general Be/XBPs. This 
indicates the neutron star might pass the outer edge of the 
equatorial disk of the Be star. 
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